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Abstract. Present day tropospheric ozone and its changes be- 
tween 1850 and 2100 are considered, analysing 15 global 
models that participated in the Atmospheric Chemistry and 
Climate Model Intercomparison Project (ACCMIP). The en- 
semble mean compares well against present day observa- 
tions. The seasonal cycle correlates well, except for some lo- 
cations in the tropical upper troposphere. Most (75 %) of the 
models are encompassed with a range of global mean tropo- 
spheric ozone column estimates from satellite data, but there 
is a suggestion of a high bias in the Northern Hemisphere 
and a low bias in the Southern Hemisphere, which could in- 
dicate deficiencies with the ozone precursor emissions. Com- 
pared to the present day ensemble mean tropospheric ozone 
burden of 337 ± 23 Tg, the ensemble mean burden for 1850 
time slice is ~30 % lower. Future changes were modelled us- 
ing emissions and climate projections from four Representa- 
tive Concentration Pathways (RCPs). Compared to 2000, the 
relative changes in the ensemble mean tropospheric ozone 
burden in 2030 (2100) for the different RCPs are: —4% 
(-16%) for RCP2.6, 2% (-7%) for RCP4.5, 1 % (-9%) 
for RCP6.0, and 7% (18%) for RCP8.5. Model agreement 
on the magnitude of the change is greatest for larger changes. 
Reductions in most precursor emissions are common across 
the RCPs and drive ozone decreases in all but RCP8.5, where 
doubled methane and a 40-150 % greater stratospheric influx 
(estimated from a subset of models) increase ozone. While 
models with a high ozone burden for the present day also 
have high ozone burdens for the other time slices, no model 
consistently predicts large or small ozone changes; i.e. the 
magnitudes of the burdens and burden changes do not appear 
to be related simply, and the models are sensitive to emis- 
sions and climate changes in different ways. Spatial patterns 
of ozone changes are well correlated across most models, but 
are notably different for models without time evolving strato- 
spheric ozone concentrations. A unified approach to ozone 
budget specifications and a rigorous investigation of the fac- 
tors that drive tropospheric ozone is recommended to help 
future studies attribute ozone changes and inter-model dif- 
ferences more clearly. 


1 Introduction 

The Atmospheric Chemistry and Climate Model Intercom- 
parison Project (ACCMIP) is designed to complement the 
climate model simulations being conducted for the Coupled 
Model Intercomparison Project (CMIP), Phase 5 (e.g. Tay- 
lor et al., 2012), and both will inform the Intergovernmental 
Panel on Climate Change (IPCC) Fifth Assessment Report 
(AR5). A primary goal of ACCMIP is to use its ensemble 
of tropospheric chemistry-climate models to investigate the 
evolution and distribution of short-lived, chemically-active 
climate forcing agents for a range of scenarios, a topic that is 
not investigated in such detail as part of CMIP5. Ozone in the 


troposphere is one such short-lived, chemically-active forc- 
ing agent, and, as it is both a pollutant and greenhouse gas, 
it straddles research communities concerned with air qual- 
ity and climate. This study is concerned with quantifying the 
evolution and distribution of tropospheric ozone in the AC- 
CMIP models, detailing the projected ozone changes since 
the pre-industrial period through to the end of the 21st cen- 
tury, with a focus on where the projected changes from the 
ensemble are robust. 

Ozone is not directly emitted and its abundance in the tro- 
posphere is determined from a balance of its budget terms: 
chemical production and influx from the stratosphere, versus 
chemical loss and deposition to the surface (e.g. Lelieveld 
and Dentener, 2000). The magnitudes of these terms are sen- 
sitive to the prevailing climate, and the levels and locations 
of ozone precursor emissions, such as nitrogen oxides (NO 
and NO 2 ; referred to as NO x ), carbon monoxide (CO) and 
volatile organic compounds (VOCs), including methane (e.g. 
Wild, 2007). A number global model studies have explored 
how changes in these drivers could affect tropospheric ozone 
abundances, from the pre-industrial period to future projec- 
tions (e.g. Johnson et al., 1999; Collins et al., 2003; Prather 
et al., 2003; Shindell et al., 2003, 2006c; Sudo et al., 2003; 
Zeng and Pyle, 2003; Mickley et al., 2004; Hauglustaine 
et al., 2005; Lamarque et al., 2005, 2011; Stevenson et al., 
2005; Brasseur et al., 2006; Dentener et al., 2006; West et 
al., 2007; Wu et al., 2008; Zeng et al., 2008; Jacobson and 
Streets, 2009; Young et al., 2009; Kawase et al., 2011). 

Estimates of past emissions of anthropogenic ozone pre- 
cursors are much lower than for the present day (e.g. Lamar- 
que et al., 2010), meaning models project large increases 
in tropospheric ozone since the pre-industrial era (Hauglus- 
taine and Brasseur, 2001; Lamarque et al., 2005; Shindell et 
al., 2006a; Cionni et al., 2011). This is in qualitative agree- 
ment with observationally based assessments (Volz and Kley, 
1988), although matching the estimated low surface ozone 
concentrations (Marenco et al., 1994; Pavelin et al., 1999) 
is challenging (Hauglustaine and Brasseur, 2001; Mickley et 
al., 2001). Future projections of anthropogenic ozone pre- 
cursor emissions used by earlier model studies often relied 
on the high population/high fossil fuel growth scenarios of 
Nakicenovic et al. (2000), meaning large emission increases 
(e.g. NO x emissions increasing nearly 4-fold between the 
1990 and 2100), resulting in large increases in tropospheric 
ozone levels (Stevenson et al., 2000; Sudo et al., 2003; Zeng 
and Pyle, 2003; Shindell et al., 2006c). However, more recent 
emission projections have included scenarios with reductions 
in anthropogenic precursor emissions (considering more ex- 
tensive air quality legislation), resulting in decreased tropo- 
spheric ozone compared to the present day (Dentener et al., 
2005; Stevenson et al., 2006; West et al., 2006, 2007, 2012; 
Kawase et al., 2011; Lamarque et al., 2011). With regard 
to natural emissions, lightning NO x emissions have gener- 
ally been thought to increase in a warmer climate (e.g. Price 
and Rind, 1994; Schumann and Huntrieser, 2007), although 
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this result is not universal (Stevenson et al., 2005; Jacob- 
son and Streets, 2009). For biogenic emissions, isoprene is 
likely the largest contributor (e.g. Guenther et al., 1995). At 
the leaf level, its emission flux depends on climate and (in- 
versely) on CO 2 concentration (Guenther et al., 2006; Ar- 
neth et al., 2010), and whether future isoprene emission is 
projected to increase (Sanderson et al., 2003; Lathieire et al., 
2005) or decrease (Arneth et al., 2007a; Young et al., 2009) 
depends on whether the CO 2 dependency is excluded or in- 
cluded (see also Pacifico et al., 2009). Biogenic emissions 
also depend on the amount and type of the vegetation, so 
projecting past and future emissions also depends on changes 
in vegetation growth, land cover and land use (Sanderson et 
al., 2003; Wiedinmyer et al., 2006; Ameth et al., 2008; Ash- 
worth et al., 2012). These changes also impact deposition 
rates (Ganzeveld et al., 2010), all together making for com- 
plex biosphere-atmosphere interactions (Arneth et al., 2010). 

The impacts of climate change on meteorology and large- 
scale atmospheric dynamics are also important for tropo- 
spheric ozone. For example, several studies report an in- 
crease in the stratospheric influx of ozone in response to 
a warming climate, resulting from a climate change-driven 
strengthening of the residual circulation (Collins et al., 2003; 
Sudo et al., 2003; SPARC-CCMVal, 2010), coupled with the 
impact of higher-than-present levels of stratospheric ozone 
(Hegglin and Shepherd, 2009; Zeng et al., 2010). On the 
other hand, increases in specific humidity in a warmer at- 
mosphere can increase the ozone loss rate, speeding up the 
reaction rate of O('D) + H 2 O (producing OH) at the expense 
of collisional quenching, 0( 1 D) + M (producing ozone again) 
(Thompson et al., 1989; Johnson et al., 1999). Higher tem- 
peratures also reduce the efficacy of peroxy acetyl nitrate 
(PAN) as a NO x reservoir, and can mean larger or smaller 
rate coefficients, depending on the activation energy of the 
given reaction. 

Despite general agreement on how the drivers impact 
global-scale shifts in tropospheric ozone, magnitudes of the 
regional changes and the overall ozone budget vary con- 
siderably between different models (Stevenson et al., 2006; 
Wild, 2007; Wu et al., 2007). With the movement to con- 
sider more physical processes and create complex Earth sys- 
tem models, uncertainty for future climate and composition 
projections may well increase (Stainforth et al., 2007). Yet, 
multi-model evaluation against present day observations and 
comparisons of the projections between models remains use- 
ful, both for benchmarking and for identifying consistent 
and contradictory results between different parameterisations 
(e.g. Dentener et al., 2006; Shindell et al., 2006b; Steven- 
son et al., 2006). In this study we analyse the multi-model 
ACCMIP ensemble ozone changes, from 1850 through to 
near- (2030) and further-term (2100) projections, using the 
latest set of scenarios developed for the CMIP5 simulations. 
This is the first study to examine the spread of modelled 
ozone responses using these scenarios, expanding the single- 
model studies of Kawase et al. (2011), Cionni et al. (2011) 


and Lamarque et al. (2011), and building on the last ma- 
jor multi-model model comparison for tropospheric ozone 
changes, coordinated by the European Union project Atmo- 
spheric Composition Change: the European Network of Ex- 
cellence (ACCENT) (Stevenson et al., 2006). The main fo- 
cus of this study is on the ensemble mean ozone change, and 
the robustness of the results across the ACCMIP ensemble. 
A detailed investigation into the “process-based” drivers of 
the ozone changes and inter-model differences is not pos- 
sible due to the limited ozone budget data and the lack of 
simulations designed to isolate particular processes (e.g. as 
per Lee et al., 2012); we recommend that this is a priority 
for future chemistry-climate model comparisons. The anal- 
ysis complements parallel investigation of the ACCMIP en- 
semble, related to climate evaluation (Lamarque et al., 2013), 
OH and methane lifetime (Naik et al., 2012; Voulgarakis et 
al., 2012), and the radiative impact of tropospheric ozone 
(Bowman et al., 2012; Shindell et al., 2012; Stevenson et al., 
2012). This study is also complementary to an investigation 
of tropospheric and stratospheric ozone in the CMIP5 models 
(Eyring et al., 2012). 

This study is organised as follows. Section 2 summarises 
the salient details of the ACCMIP models and the simula- 
tions analysed, followed by a comparison of the model emis- 
sions in Sect. 3. Section 4 discusses the present day distribu- 
tion of tropospheric ozone and the inter-model differences, 
and presents a reference comparison of the models against a 
range of ozonesonde and satellite -based measurement data. 
The modelled ozone changes for the different scenarios are 
documented in Sect. 5, followed by a brief discussion of all 
the results in Sect. 6. Finally, Sect. 7 summarises the main 
conclusions and recommendations for future multi-model in- 
vestigations of tropospheric ozone. 


2 Models, simulations and analysis details 

Here we provide brief details of the ACCMIP models and 
simulations, together with some details on the analysis per- 
formed in this study. Lamarque et al. (2013) provide a more 
complete description of the models, including appropriate 
references, and further details for the simulations. 

2.1 ACCMIP models 

Table 1 summarises the models, scenarios and their time pe- 
riods analysed in this study (the tropospheric ozone burdens 
are discussed in later sections of the text). For this study, we 
used the output from 15 models, although not all of them 
provided output for every scenario and period, as indicated 
by ” in Table 1. Note that the analysis here does not in- 
clude the ACCMIP model NCAR-CAM5.1 as this did not 
calculate ozone. 

Most of the ACCMIP models are climate models with 
atmospheric chemistry modules, run in atmosphere-only 
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Table 1 . Tropospheric ozone burdens (Tg) for the individual models and the ACCMIP mean. Also shows which simulations and time slices 
are available from each model for this study. Not all relevant variables are available for each model/time slice. 


Model 


Hist 


RCP2.6 

RCP4.5 

RCP6.0 

RCP8.5 


1850 

1980 

2000 

2030 

2100 

2030 

2100 

2030 

2100 

2030 

2100 

CESM-CAM-superfast 

192 

288 

302 

278 

230 

- 

- 

289 

252 

328 

384 

CICERO-OsloCTM2 a 

206 

287 

308 

296 

247 

312 

274 

- 

- 

326 

343 

CMAM 

239 

310 

323 

307 

266 

330 

293 

- 

- 

342 

371 

EMAC 

259 

352 

378 

- 

- 

379 

342 

- 

- 

399 

441 

GEOSCCM 

250 

333 

346 

- 

- 

- 

- 

- 

- 

- 

- 

GFDL-AM3 

264 

370 

378 

367 

317 

389 

356 

390 

359 

410 

484 

GISS-E2-R b 

252 c 

337 

344 

341 

304 e 

352 

32 l e 

352 

339 e 

379 

417 e 

GISS-E2-R-TOMAS 

261 

350 

359 

- 

- 

- 

- 

- 

- 

- 

- 

HadGEM2 

227 

289 

307 

303 

262 

316 

295 

- 

- 

330 

377 

LMDz-OR-INCA b 

247 c 

322 

339 d 

321 

278 e 

342 

310 e 

329 

306 e 

351 

374 e 

MIROC-CHEM 

239 

321 

341 

325 

283 

- 

- 

338 

304 

356 

374 

MOCAGE 

272 

322 

327 

323 

299 

- 

- 

333 

358 

400 


NCAR-CAM3.5 

221 

318 

336 

317 

263 

336 

294 

322 

285 

349 

386 

STOC-HadAM3 

234 

332 

348 

329 

272 

- 

- 

- 

- 

367 

385 

UM-CAM 

226 

304 

322 

323 

293 

338 

322 

- 

- 

351 

397 

Mean 

239 

322 

337 

319 

276 

344 

312 

336 

309 

357 

395 

Sdev ( % of mean) 

22 (9%) 

24 (8%) 

23 (7%) 

22 (7%) 

25 (9%) 

26 (8%) 

26(8%) 

31 (9%) 

35 (11%) 

26 (7%) 

36(9%) 


a Simulations are a single year. 

b These models submitted transient simulations. Their “time slice” means represent 10-yr averages about the given decade, except as noted. 
c Mean of 1850-1 859. 

d Mean of 1996-2000 (when the transient simulation stops). 
e Mean of 2091-2100. 


mode; i.e. the models are driven by sea-surface temper- 
ature (SST) and sea-ice concentrations (SICs). GISS-E2- 
R uniquely was run as a fully coupled ocean-atmosphere 
climate model, although the closely related GISS-E2-R- 
TOMAS model was run with SSTs and SICs prescribed. 
CICERO-OsloCTM2 and MOCAGE are chemical transport 
models (CTMs), with MOCAGE using offline meteorologi- 
cal fields from an appropriate simulation of a climate model, 
and CICERO-OsloCTM2 using offline meteorological fields 
from a single year of a reanalysis dataset. Except for the 
CTMs, and LMDz-OR-INCA, STOC-HadAM3 and UM- 
CAM, the calculated ozone concentrations are used in the 
climate model radiation code, making most of the models 
chemistry-climate models (CCMs). 

The model chemical schemes vary greatly in their com- 
plexity (e.g. as measured by the number of species and re- 
actions), particularly in the range of non-methane VOCs 
(NMVOCs) that they simulate. Complexity ranges from 
the simplified and parameterized schemes of CMAM (no 
NMVOCs) and CESM-CAM-superfast (isoprene as the only 
NMVOC), to the intermediate schemes of HadGEM2 and 
UM-CAM (include < C 3 -alkanes), to the more complex 
schemes of the other models, which include the more reac- 
tive, chiefly anthropogenic NMVOCs (e.g. higher alkanes, 
alkenes and aromatic species), as well as lumped monoter- 
penes. Some representation of stratospheric chemistry is in- 
cluded in many models, with the exception of CICERO- 
OsloCTM2, HadGEM2, LMDz-OR-INCA, STOC-HadAM3 
and UM-CAM. CICERO-OsloCTM2 uses monthly-varying 


climatological ozone values from a previous model sim- 
ulation (except for the lowest ~2.5 km of the strato- 
sphere), LMDz-OR-INCA uses a constant (in time) strato- 
spheric ozone climatology (Li and Shine, 1995), whereas the 
other models without detailed stratospheric chemistry used 
the time varying stratospheric ozone dataset of Cionni et 
al. (2011). 

2.2 Scenarios and time slices 

The ACCMIP simulations broadly correspond with the 
CMIP5 scenarios (Taylor et al., 2012). Historical (hereafter 
Hist) simulations cover the preindustrial period to the present 
day, while a range of Representative Concentration Pathways 
(RCPs) (van Vuuren et al., 2011) cover 21st century pro- 
jections. These latter scenarios are named for their nominal 
radiative forcing level (2100 compared to 1750), such that 
RCP2.6 corresponds to 2.6 Wm" 2 , RCP4.5 to 4.5 Wm" 2 , 
RCP6.0 to 6.0 Wm“ 2 and RCP8.5 to 8.5 Wm“ 2 . Ozone pre- 
cursor emissions from anthropogenic and biomass burning 
sources were taken from those compiled by Lamarque et 
al. (2010) for the Hist simulations, whereas emissions for 
the RCP simulations are described by Lamarque et al. (2013) 
(see also Lamarque et al., 2011; van Vuuren et al., 2011). 
Excluding methane emissions, all the RCPs include reduc- 
tions and redistributions of ozone precursor emissions mov- 
ing through the 21st century. Natural emissions, such as CO 
and VOCs from vegetation and oceans, and NO x from soils 
and lighting, were determined by each model group. The 
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emissions used by the individual models are discussed fur- 
ther in Sect. 3. 

With the exception of GISS-E2-R and LMDz-OR-INCA, 
each model conducted a set of time slice simulations for each 
scenario. In this study, we analyse output from the 1850, 
1980 and 2000 time slices from the Hist scenario, and the 
2030 and 2100 time slices for the RCPs, to provide near- 
term and longer-term perspectives. Except for the CTMs 
and GISS-E2-R, each model used climatological SSTs and 
SICs from coupled ocean-atmosphere CMIP5 simulations of 
a closely related climate model, typically averaged for the 
10 yr about each time slice (e.g. 2026-2035 for the 2030 
time slice), although some models had interannually varying 
boundary conditions. CICERO-OsloCTM2 used the same 
meteorology for each simulation, whereas MOCAGE was 
run with meteorological fields from a climate model running 
the appropriate time slice and scenario, rather than directly 
using the SSTs and sea-ice to drive an atmosphere model. 
The number of years that the ACCMIP models simulated 
for each time slice mostly varied between 4 and 12 yr for 
each model, although CICERO-OsloCTM2 only simulated a 
single year. However, as the boundary conditions (including 
biomass burning emissions) were constant for each year of a 
given time slice for most models, “interannual” variability is 
generally small (see Sect. 4). 

GISS-E2-R and LMDz-OR-INCA both conducted tran- 
sient simulations, and the data analysed in this study were 
averaged for the decade about each time slice (e.g. 1 976— 
1985 for the 1980 time slice), with some minor exceptions as 
noted in Table 1 . 

2.3 Analyses: tropopause definition and statistical 
definitions 

Throughout this analysis, the troposphere is defined as air 
with ozone concentrations less than or equal to 150ppbv 
(Prather et al., 2001), which is simple to employ and allows 
comparison against other model studies (e.g. Stevenson et 
al., 2006). For a consistent definition of the troposphere for 
all time slices, the definition is applied using the ozone from 
the Hist 1850 time slice mean, applied on a per model basis, 
and varying by month. We used the 1850 time slice to avoid 
issues with different degrees of stratospheric ozone depletion 
across the ensemble, particularly in the Southern Hemisphere 
(SH). While fixing the definition means we compare a consis- 
tent region of the atmosphere between different time slices, 
it does ignore the fact the tropopause height will likely alter 
with climate change (Santer et al., 2003a, b). Furthermore, 
values for the tropospheric ozone burdens and columns are 
obviously sensitive to the tropopause definition (Wild, 2007; 
Prather et al., 2011), and the differences generally amount to 
± 5 % compared to using the 2000 time slice mean ozone, al- 
though they can exceed 10 % for a few models for the RCPs. 
We apply a pressure based tropopause definition to ensure a 


consistent comparison of the modelled tropospheric column 
ozone against satellite data in Sect. 4.2. 

A statistical analysis of whether the projected ozone 
changes are significant against interannual variability is not 
possible with the ACCMIP data, mainly because the inter- 
annual variability is insufficiently characterised by the time 
slice runs, as most have constant SSTs and SICs on a year- 
to-year basis, and all have constant biomass burning emis- 
sions. Instead, we assess whether a given multi-model mean 
ozone change is significantly different from zero by using a 
paired sample Student’s t test (e.g. Wilks, 2006). Changes are 
considered significant if the (absolute) mean change from all 
the models is greater than 2 times the standard error for the 
mean change (i.e. approximately the 5 % level). One weak- 
ness of this analysis is that it cannot highlight regions where 
the models agree that the changes are not significant with re- 
spect to interannual variability (see Tebaldi et al., 2011). 

For the most part, output from the models was interpo- 
lated to the grid used by Cionni et al. (2011), who compiled 
the ozone dataset recommended for use in the CMIP5 sim- 
ulations (5° by 5° latitude/longitude and 24 pressure levels). 
However, the tropospheric burden and columns were calcu- 
lated on a model’s native grid. 

3 Emissions: differences and similarities between 
models 

While one goal of ACCMIP was for models to match each 
other’s ozone precursor inputs as closely as possible, differ- 
ences in model parameterisations and complexity means that 
some model diversity is unavoidable. In particular, natural 
emissions were not prescribed as part of the experimental de- 
sign, and their differing treatment between models broadens 
the range of ozone precursors. Such differences are examples 
of why we need model comparisons. 

Figure 1 shows the range in ozone precursor emissions in 
the ACCMIP models, presenting box-whisker plots for each 
scenario and time slice. The number of models that constitute 
the spread of the data is different for different simulations - 
see Table 1 . Tabulated emission data for the individual mod- 
els can be found in Table S 1 in the Supplement. 

Figure la shows that the tropospheric methane burden is 
generally well constrained in the simulations, with the in- 
terquartile range (IQR) being 3-5.5% of the mean burden. 
The close agreement is due to all models except LMDz- 
OR-INCA having used prescribed methane surface con- 
centrations for the Hist simulations, and only GISS-E2-R 
and LMDz-OR-INCA not prescribing concentrations for the 
RCP simulations. The methane burden approximately dou- 
bles from 1850 to 2000, but 2100 burdens are 30 %, 10 % and 
2.5 % lower than 2000, for the RCP2.6, RCP4.5 and RCP6.0 
scenarios respectively. For RCP8.5, by 2100 the methane 
burden has more than doubled again compared to 2000. 
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(a) Methane burden 


(e) Total VOC 



(b) CO 


(c) Total NOx 


1 850 1 980 2000 2030 21 00 2030 21 00 2030 21 00 2030 21 00 
time slice 



time slice 

(f) Biogenic VOC 
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Fig. 1 . Burdens and emissions of ozone precursors used in the ACCMIP scenarios and time slices, showing (a) the tropospheric methane 
burden, and yearly total (b) CO emissions, (c) total NO x emissions, (d) lightning NO x emissions and (e) total VOC emissions, which are 
then split into (f) biogenic VOCs and (g) other VOCs. The spread of the emissions/burden in each model is indicated for each scenario and 
time slice, with filled box showing the interquartile range, the whiskers indicating the full range, and the dot and line indicating the mean and 
median respectively. Not all models completed all the scenarios: the number of models used to determine each box/whisker is indicated. 


The general trends in CO (Fig. lb) and total NO x (Fig. lc) 
emissions are similar, increasing from 1850 to 2000 for the 
Hist scenario, decreasing thereafter for all the RCPs, ex- 
cept for the 5% higher NO x emissions for the 2030 time 
slice of RCP8.5, compared to Hist 2000. NO x emissions 
show a stronger increase over the 20th century, with the 
mean trebling between 1850 and 2000, whereas CO emis- 


sions slightly more than double. Across all RCPs, by 2100 
CO and NO x emissions are lower by 30-45 % compared to 
2000 . 

Both CO and NO x emissions show a greater degree of vari- 
ation between the models than the methane burden. For a 
given time slice, the IQRs vary between approximately 10- 
30% of the corresponding mean emission, whereas the full 
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range (maximum minus minimum emission) is between 20- 
100% of the mean emission. The spread is due to the vary- 
ing natural emissions (NO x from soils and lightning; CO 
from oceans and vegetation), as well as the less complex 
chemistry schemes used in some models including extra CO 
emissions as a surrogate for missing NMVOCs (e.g. CMAM, 
HadGEM2). 

An example of the variation in natural emissions can be 
seen from Fig. Id, which shows the spread in the lightning 
NO x source (LNO x ) for the models. Parameterisations for 
LNO x are generally dependent on cloud top heights and con- 
vective mass fluxes (e.g. Price and Rind, 1992; Allen and 
Pickering, 2002), which likely show large variability be- 
tween models, accounting for spread. The IQRs are gen- 
erally 40-55 % of the mean emission, and the full range 
is 90-170% of the mean. The maximum emissions come 
from the MIROC-CHEM, whose LNO x was set mistakenly 
high (by 60%) in the ACCMIP simulations. The minimum 
emissions (< 2 TgNyr -1 ) are generally from the HadGEM2 
LNO x (who also mistakenly implemented LNO x ), although 
emissions are also low for CMAM for the RCP simula- 
tions. Our knowledge of LNO x is generally poor, but (ex- 
cluding HadGEM2 and MIROC-CHEM) LNO x for the Hist 
2000 simulation ranges between 3. 8-7. 7 TgNyr -1 , within 
the range of 5 ±3 TgNyr -1 estimated by Schumann and 
Huntrieser (2007) for a range of LNO x parameterisations. 

Possible changes in lightning activity with climate change 
(e.g. Williams, 2009 and refs, therein) have been recognised 
as potentially important for LNO x and the subsequent im- 
pact on tropospheric composition (e.g. Price and Rind, 1994; 
Hauglustaine et al., 2005; Fiore et al., 2006; Schumann and 
Huntrieser, 2007; Zeng et al., 2008), even if only the spatial 
distribution changes (Stevenson et al., 2005). An increase in 
LNO x from 2000 to 2100 (RCP8.5; strongest warming) is 
generally robust across the ACCMIP models, and ranges in 
magnitude from 10-75 %. CMAM is an outlier in this case, 
with 45% lower emissions for RCP8.5 2100 compared to 
Hist 2000. CMAM is also the only model using an LNO x 
parameterisation based on the study of Allen and Picker- 
ing (2002). Jacobson and Streets (2009) also modelled lower 
LNO x in a warmer climate, using a different parameterisa- 
tion based on cloud ice. Clearly further study is required into 
the implications of the use of different parameterisations for 
LNO x , and the different sensitivities across models. 

Figure le shows the total VOC emissions used in the AC- 
CMIP models. As with LNO x , the emissions cover a wide 
range and there is no clear trend in the mean across the 
simulations. Many of the reasons for the differences are fa- 
miliar from the above discussion, particularly the fact that 
some models include more VOC species than others. An- 
other reason for the spread in VOC emissions comes from 
treatment of VOCs of biogenic origin, particularly isoprene, 
which likely dominates the total NMVOC emissions (e.g. 
Guenther et al., 1995). Figure If shows the spread of bio- 
genic VOC emissions between the ACCMIP models; these 


emissions are mostly isoprene. EMAC, GEOSCCM, GISS- 
E2-R and STOC-HadAM3 simulations were the only ones to 
include climate-sensitive isoprene emissions, and these are 
the only models with a positive change in VOC emissions 
between Hist 2000 and RCP8.5 2100, arising from the posi- 
tive temperature dependence of isoprene emission algorithms 
(e.g. Guenther et al., 2006; Arneth et al., 2007b). Arneth et 
al. (2011) noted that the isoprene emission computed from a 
given algorithm is sensitive to the input meteorological data 
and vegetation boundary conditions, giving further cause for 
variation in VOC emissions between models. For the rest 
of the ACCMIP ensemble, if they included isoprene chem- 
istry, constant present day isoprene emissions were used for 
all simulations. Figure lg shows that the trend in the non- 
biogenic VOC emissions (anthropogenic plus biomass burn- 
ing) broadly resembles that of CO and NO x , albeit with the 
range of VOC complexity resulting in a larger spread of the 
emission total between the models (IQR is 30-100% of the 
mean emission). 


4 Present day ozone distribution and 
model-observation comparison 

This section presents the tropospheric ozone burden, the tro- 
pospheric ozone budget (for a limited subset of models), and 
the distribution of tropospheric ozone (surface, column and 
zonal mean) as simulated by the ACCMIP models for the 
Hist 2000 simulation. The ozone concentrations and columns 
from this simulation are then compared against observational 
datasets, from both satellites and ozonesondes. We do not 
over-interpret these comparisons since they exclude observa- 
tional error, and they differ in the periods covered. The em- 
phasis is largely on the distribution and measurement-model 
comparison for the ensemble mean, describing the spread 
of model results with statistical metrics. The distributions of 
ozone for all the individual models can be found in the sup- 
plementary material (Figs. S1-S3). 

4.1 Tropospheric ozone burden and budget 

Figure 2a shows the annual mean tropospheric ozone bur- 
den for all the models for Hist 2000, as well as the AC- 
CMIP mean. Values for the tropospheric burden for all mod- 
els and scenarios can also be found in Table 1. The mean 
burden is 337±23Tg, very close to the 336±27Tg re- 
ported for a subset of the ACCENT models (Stevenson et al., 
2006) (344 ± 39 Tg for the full ACCENT ensemble), and the 
335±10Tg estimated from measurement climatologies by 
Wild (2007), although the latter estimate is from pre-2000 
ozone data. Figure 2a also indicates the uncertainty in the 
ozone burden, as represented by the standard deviation of 
the range of burdens computed for individual years of the 
time slice. As might be expected from successive years of the 
same boundary conditions (for most models), the uncertainty 
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Table 2. Tropospheric ozone budget statistics for a subset of ACCMIP models for the Hist 2000 time slice, showing chemical production 
(P), chemical loss (L), deposition (D), the inferred stratospheric influx (S), and the lifetime (r). 


Model 


Flux terms/Tg 

ozone a 1 


r/days 


P 

L 

D 

S 


CESM-CAM-superfast 

3877 

3638 

687 

448 

25.5 

GEOSCCM 

4692 

3853* 

1240 

401 

24.8 

GFDL-AM3 

5853 

5089 

1240 

476 

21.8 

NCAR-CAM3.5 

4494 

4112 

842 

460 

24.8 

STOC-HadAM3 

5989 

5050 

1350 

411 

19.9 

UM-CAM 

4358 

3816 

1205 

663 

23.4 

ACCENT mean (± sdev) 

5 110 ±606 

4668 ±727 

1003 ± 200 

552 ± 168 

22.3 ±2.0 


* Loss flux includes wet deposition of oxidised nitrogen compounds. 




Fig. 2. (a) Tropospheric ozone burdens for the ACCMIP models 
from the Hist 2000 simulations. Error bars for the models indicate 
the variability in the burden between different years of a model’s 
time slice (± 1 std. dev.). The error bar on the ensemble mean bur- 
den indicates the inter-model spread of the burden (± 1 std. dev.), 
(b) Distribution of the mean ozone burden throughout the tropo- 
sphere for the mean model, using the “boxes” recommended by 
Lawrence et al. (2001) for reporting OH concentrations. The solid 
line indicates the tropopause (see text for definition). 


is small, and the standard deviations are less than 2 % of the 
burden. 

There is a significant correlation (a- = 0.67) between the 
modelled ozone burden and the total VOC emissions. With 
the spread of VOC emission and treatment between the mod- 
els, it is difficult to rationalise this correlation satisfactorily, 
although Wild (2007) demonstrated increased VOC emis- 
sions lead to an increased ozone burden. There is not a simi- 
lar correlation between the ozone burden and total NO x emis- 
sions. 

Figure 2b indicates the distribution of the mean ozone bur- 
den throughout the troposphere, using the regions defined by 
Lawrence et al. (2001) (to describe the distribution of OFI) 
to give a mass-weighted view of the zonal ozone distribu- 
tion. The hemispheric asymmetry in ozone is apparent from 
Fig. 2b, which shows that 57.5% of the ozone mass is in 
the NFI. The NFI extratropics has 60 % more ozone than the 
SFI extratropics overall, but the NFI tropics has only slightly 
more ozone (~3 %) than the SFI tropics overall. The greatest 
burdens are found in the extratropical upper troposphere, re- 
flecting the importance of stratosphere-to-troposphere trans- 
port of ozone. The next largest ozone burdens are found in the 
comparatively more polluted NFI lower troposphere, as well 
as the tropical upper troposphere. This latter region is im- 
pacted by convective transport of ozone precursors and light- 
ning emissions (Jacob et al., 1996). The standard deviation in 
the fractional distribution of ozone is also in Fig. 2b, show- 
ing that the model uncertainty in distribution of ozone mass 
is largest in the NFI extratropics and in the upper troposphere 
in general, consistent with the results described in Sect. 4.2. 

Table 2 presents the present day tropospheric ozone bud- 
get terms for the six models for which there are sufficient 
data. We do not report the ensemble mean result due to both 
the limited number of models and that the chemical produc- 
tion (P) and loss (L) terms were not calculated in a consis- 
tent manner (e.g. whether ozone loss through oxidised ni- 
trogen species was considered); ozone dry deposition (D) 
was calculated consistently for these models however. Fol- 
lowing Stevenson et al. (2006), the net influx of ozone from 
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Fig. 3. ACCMIP ensemble mean, annual mean ozone climatologies and their inter-model variability, for the 2000 time slice of the Historical 
simulation. Top row shows zonal mean ozone, the middle row shows the tropospheric ozone column, and the bottom row shows surface 
ozone. For each row, the left hand panels show the absolute values of the ozone variable: ppbv for the zonal mean and surface concentrations, 
and Dobson units (DU) for the tropospheric column. The middle panels show the absolute standard deviations in the same units. The right 
hand column expresses the standard deviation as a percentage of the ensemble mean value (also known as the coefficient of variation). Note 
that each panel has a different scale. 


the stratosphere (S; n f) is calculated from the residual of the 
other terns (Si n f=P - L + D), and the tropospheric ozone 
lifetime (r) is calculated using the burdens (B) in Table 1 
(r = B/F, where F = L + D = P + Sj n f). As with the burden, all 
flux terms are defined using the 150 ppbv ozone contour as 
the tropopause (from Hist 1850). 

Several differences are apparent from comparing the AC- 
CMIP results against the ACCENT study (the ACCENT 
mean data are shown in Table 2), although we caution that, 
with the limited amount of ACCMIP data, generalisations 
about how the modelled budget has changed since ACCENT 
are hard to make. For GFDL-AM3 and STOC-HadAM3, P 
is much higher than the ACCENT mean, whereas P is much 
lower for CESM-CAM-superfast, NCAR-CAM3.5 andUM- 
CAM. For L, the ACCMIP models are broadly ordered the 
same as P, although - unlike P - all the L terms all sit 
within the range described by the ACCENT mean and stan- 
dard deviation. The ACCMIP models with the lower P and 
L have lower total VOC emissions (see Table Sib), which 
could go some way to explaining the range in Table 2 (e.g. 
Wild, 2007). These models also have the longest tropospheric 
ozone lifetimes. For D, the ACCMIP results span nearly a 


factor of two between CESM-CAM-superfast and STOC- 
HadAM3. The fact that D does not simply correlate with 
B underlines that there are differences in the ozone spatial 
distribution and the deposition implementation between the 
models (see Lamarque et al., 2013), which has implications 
for assessing the impacts of ozone concentration changes on 
the biosphere (e.g. Sitch et al., 2007). For S; n f, all the AC- 
CMIP models are encompassed by the ACCENT mean and 
standard deviation, and, furthermore, S; n f for the six models 
is within the 540 ± 140 Tg yr _ 1 range suggested using obser- 
vational constraints of stratosphere-to-troposhere exchange 
(Olsen et al., 2001; Wild, 2007). However, determining the 
net stratospheric influx by budget closure will likely give a 
different value to that calculated using transport diagnostics 
within the model (e.g. Sanderson et al., 2007), and - as with 
all the budget terms - there will be some sensitivity to the 
tropopause definition (Wild, 2007). 


www.atmos-chem-phys.net/13/2063/2013/ 


Atmos. Chem. Phys., 13, 2063-2090, 2013 






2072 


P. J. Young et al.: Pre-industrial to end 21st century projections of tropospheric ozone 


4.2 Zonal mean, tropospheric column and surface 
ozone from Hist 2000 

Figure 3 shows the ensemble mean, annual mean distribu- 
tion of ozone, presenting the zonal mean, tropospheric col- 
umn and surface ozone concentrations, as well as their inter- 
model variability. The general patterns of the ozone distribu- 
tion in Figs. 3a, 3d and 3g are consistent with those reported 
from satellite (Fishman et al., 1990; Ziemke et al., 2011) and 
ozonesonde (Logan, 1999; Thompson et al., 2003) measure- 
ments, as well as the multi-model data shown by Stevenson et 
al. (2006). The increase in ozone concentration with height is 
clear, in accordance with the increase in ozone lifetime. Con- 
vective lifting of low-ozone air masses coupled with lofting 
of ozone precursors (Lawrence et al., 2003; Doherty et al., 
2005) results in the characteristic tropical zonal mean profile 
in Fig. 3a. The hemispheric asymmetry in mid-tropospheric 
ozone concentrations reflects the larger input of stratospheric 
ozone in the NH, due to the stronger Brewer-Dobson circula- 
tion there (Rosenlof, 1995), as well as the larger emissions 
of ozone precursors (Lamarque et al., 2010). While both 
the tropospheric column (Fig. 3d) and surface concentrations 
(Fig. 3g) also show higher ozone levels over ozone precur- 
sor source regions, the plots also indicate enhanced concen- 
trations downwind of the these regions, due to transport of 
ozone and ozone precursors, including “reservoir” species, 
such as PAN (Moxim et al., 1996; Fiore et al., 2009). Fig- 
ure 3d also shows the “wave-1” pattern in the tropical tro- 
pospheric ozone column (Thompson et al., 2003; Ziemke et 
al., 2011), with a minimum in ozone over the Pacific Ocean 
and maximum over the Atlantic. Surface ozone concentra- 
tions are also very low over the equatorial Pacific Ocean. 

There is generally good agreement between the models for 
the zonal mean profile of ozone. Figure 3c shows that the 
standard deviation is less than 20 % of the mean throughout 
much of the troposphere, with exception of some lower tro- 
posphere regions and the upper troposphere. The spatial pat- 
terns of the spread in surface ozone concentrations in Fig. 3h 
and i suggest that much of the lower troposphere variability is 
over regions with large anthropogenic, pyrogenic or biogenic 
emissions, where both the absolute and relative uncertainty 
is largest. In anthropogenic and biomass burning source re- 
gions, much of the model diversity could reflect the spread 
in VOC composition (low vs. high reactivity species), which 
means different ozone production efficiencies (e.g. Russell 
et al., 1995), as well as using different injection heights for 
biomass burning emissions. For tropical Africa and espe- 
cially South America, large variations are apparent over iso- 
prene source regions (surface and column), which reflects 
differences in the total emission (some models have no iso- 
prene), as well as potentially differences in isoprene chem- 
istry (Archibald et al., 2010). 

Large model variation is also found for the high latitude 
SH, chiefly for the tropospheric column. Tropospheric ozone 
levels in the SH are generally low, but there is relatively 


large diversity in the overhead stratospheric ozone column 
(standard deviation for the total ozone column is 10-15 % of 
the mean; not shown). This results in a spread in the strato- 
spheric input as well as potentially some impacts through 
changes in photolysis rates, for those models with photolysis 
schemes that use the model-calculated ozone column (e.g. 
Fuglestvedt et al., 1994; see also Voulgarakis et al., 2012). 
There is less (relative) variation in the tropospheric column 
in the Northern Hemisphere (NH), coupled with less spread 
between models for the total ozone column (not shown). But 
larger uncertainty for the surface at high latitudes could be re- 
lated to differences in precursor transport and chemistry from 
lower latitudes (Eckhardt et al., 2003; Shindell et al., 2008; 
Christoudias et al., 2012). As was also found by Stevenson 
et al. (2006), there are large relative uncertainties in tropo- 
spheric column ozone over the equatorial Pacific Ocean, but 
the concentrations here are very low. 

4.3 Comparison to ozonesondes and satellite data 

Figure 4 compares the ACCMIP mean, median and indi- 
vidually modelled ozone concentrations from the Hist 2000 
simulation against ozonesonde data, in the same manner 
as Stevenson et al. (2006). Ozonesonde measurements are 
taken from datasets described by Logan (1999) (representa- 
tive of 1980-1993) and Thompson et al. (2003) (represen- 
tative of 1997-2011), and consist of 48 stations, split 5, 15, 
10 and 18 between the SH extratropics, SH tropics, NH trop- 
ics and NH extratropics respectively. The models were sam- 
pled at the ozonesonde locations. In addition, Fig. 4 shows 
satellite-derived ozone concentrations retrieved from the Tro- 
pospheric Emission Spectrometer (TES), which is a Fourier 
transform spectrometer on board the NASA Aura spacecraft 
in 2004 (Beer, 2006). TES ozone from a 2005-2010 climatol- 
ogy was interpolated to the same grid as the ACCMIP mod- 
els and sampled at the ozonesonde locations. Figure 4 also 
shows the ACCENT model mean (Stevenson et al., 2006), to 
place the ACCMIP results in context of recent multi-model 
comparisons. The correlation and mean normalised bias er- 
ror (MNBE) are shown for multi-model mean from the AC- 
CMIP and ACCENT ensembles, relative to the ozonesonde 
observations. 

Both the ACCMIP ensemble mean and median are within 
the standard deviation of the observations for most loca- 
tions and altitudes, with the winter NH extratropical com- 
parison being a notable exception. Indeed, compared to 
the mean observations, the largest relative errors are found 
for the NH extratropics, where the mean model overesti- 
mates the concentrations, and SH tropics, where the mean 
model underestimates the concentrations. The individual 
model biases in these locations are significantly correlated 
with total VOC emissions (r = 0.57; i.e. more VOC emis- 
sions give a more positive, or less negative, bias), although 
several other chemical and transport factors likely play a 
role. However, the mean model captures the annual cycle 
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Fig. 4. Comparison of the annual cycle of ozone, between ozonesonde observations (black circles) and the ACCMIP ensemble mean (solid 
red line), ACCMIP ensemble median (dashed red line) and the ACCENT ensemble mean (blue line) (Stevenson et al., 2006). Ozone concen- 
trations from TES (purple line) are also shown. ACCMIP model data is from the 2000 time slice mean of the Historical experiment. Model 
and observational data were grouped into four latitude bands (90° S to 30° S, 30° S to 0°, 0° to 30° N and 30° N to 90° N) and sampled at 
three altitudes (700 hPa, 500 hPa and 250 hPa), with the models sampled at the ozonesonde locations before averaging together. The individ- 
ual ACCMIP models are represented by the thin grey lines, with the grey shaded area indicating ± 1 standard deviation about the ACCMIP 
ensemble mean, showing the model spread. Error bars on the observations indicate the average interannual standard deviation for each group 
of stations. The correlation (r) and mean normalised bias error (mnbe) for the ACCMIP (red) and ACCENT (blue) ensemble means versus 
the observations are also indicated in each panel. This figure is an update of Fig. 2 of Stevenson et al. (2006). 


in ozone concentrations extremely well in most locations 
(as measured by the correlation coefficient), suggesting that, 
broadly speaking, the seasonality in circulation patterns, 
stratosphere-to-troposphere exchange and natural emissions 
(chiefly biomass burning in the tropics, and isoprene in the 
NH extratropics) is captured well. The statistics for the NH 
tropical mid and upper troposphere suggest that the season- 
ality is less well modelled, although we note that, (1) the 
observed-model correlation is significant (r > 0.58), (2) there 
is considerable observed interannual variability in ozone the 
upper troposphere, and (3) the bias and correlation are im- 
proved compared to the ACCENT mean. Compared to AC- 
CENT, the correlation is improved with the ACCMIP mean 
model for most locations, and the bias for some locations, 
both most prominently in the NH. 

Except for the NH Tropics at 250 hPa, the TES data are 
positively biased compared to the ozonesondes, broadly con- 


sistent with the 2-10 ppbv high bias that Nasser et al. (2008) 
noted for TES (see also Zhang et al., 2010). However, taking 
the interannual variability into account, and the fact that we 
have neither applied the TES operators in this analysis (Wor- 
den et al., 2007), nor considered measurement uncertainty, 
the TES and ozonesonde data are not notably different. This 
positive bias means that, compared to the ozonesonde data, 
the ACCMIP mean model bias against TES is improved for 
the NH extratropics, about the same for the NH tropics (op- 
posite in sign), but worsened for the SH; changes in cor- 
relation compared to the ACCMIP-ozonesonde comparison 
are marginal. As mentioned above both the ozonesondes and 
TES see a sharp increase in ozone between March and April 
at EQ-30° N, not captured in the ACCMIP mean. TES (and 
several ACCMIP models) do not show the same low values 
in the winter months as the ozonesonde data. For TES, this 
is likely due to lower thermal contrast which will make the 
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Fig. 5. Normalised biases for the ACCMIP models (Hist 2000 simu- 
lation) against the ozonesonde measurements compiled into regions 
by Tilmes et al. (2012). Each region is colour-coded, and the con- 
stituent ozonesonde sites are indicated in the top panel. Biases are 
shown for (bottom to top) 750, 500 and 250 hPa. The box, whiskers, 
line and dot show the interquartile range, full range, median and 
mean biases respectively, in a similar style to Fig. 1. The y-axis has 
the same scale in each panel. 


satellite retrievals relax to an a priori value (Bowman et al., 
2006). Bowman et al. (2012) pursue comparisons of the AC- 
CMIP models against TES further. 

Figure 5 makes a similar comparison to ozonesonde data, 
this time using the compilation of Tilmes et al. (2012). This 
dataset mostly consists of the same station data described 
by Logan (1999) and Thompson et al. (2003), but covering 
1995-2009, and aggregated into 12 regions that exhibit sim- 
ilar ozone concentration characteristics (see the top panel of 
Fig. 5 and Tilmes et al., 2012). The figure presents the mean, 
median and spread of the MNBE for the individual ACCMIP 
models (cf. Fig. 1), showing that the full range of perfor- 
mance encompasses positive and negative biases for each re- 
gion and altitude. 

The information in Fig. 5 is consistent with that in Fig. 4, 
but with more longitudinal information. For instance, we see 
that the negative bias in the SH tropical ozone is driven 
by the less favourable comparison of the model mean with 
the sites in the Atlantic/ Africa region (dark green), and the 
sign of the bias is consistent across more than 75% of the 
models. A positive bias is apparent in all the NH extratrop- 
ical regions in the low and mid-troposphere, and again is 
shared by the majority of the models. Figure 5 also shows 
low biases for the high latitude regions in the upper tro- 
posphere/lower stratosphere (UT/LS; a region not shown in 
Fig. 4). A comparison of the ACCMIP mean total ozone col- 
umn against satellite measurements from the merged Total 
Ozone Mapping Spectrometer/solar backscatter ultraviolet 
(TOMS/SBUV) data (Stolarski and Frith, 2006) suggests that 
the models overestimate the total ozone column by around 
5 % at high latitudes (not shown), although the total column 
bias is not necessarily related to UT/LS biases. Validation of 
stratospheric ozone in these models is beyond the scope of 
this study, but this would help resolve whether ozonesonde- 
model comparisons at higher altitudes are consistent with the 
satellite data. (A comparison of the ensemble mean ozone 
column against TOMS data can be found in the supplemen- 
tary material, along with a comparison of late twentieth cen- 
tury trends.) 

Tropospheric ozone columns are available from a combi- 
nation of the Ozone Monitoring Instrument (OMI) and Mi- 
crowave Limb Sounder (MLS) data. Figure 6 compares the 
ACCMIP mean tropospheric column ozone against the OMI- 
MLS climatology derived by Ziemke et al. (2011), covering 
October 2004 to December 2010. Tables 3 and 4 summarise 
the comparisons between the OMI/MLS data and individual 
models, showing the global (60° S-60° N) column biases and 
spatial correlations, and biases by latitude bands respectively. 
For each model, the column was defined using the tropopause 
pressures provided by Ziemke et al. (2011) (from the Na- 
tional Centers for Environmental Prediction, NCEP), mean- 
ing that Fig. 6a is subtly different from Fig. 2d. 

From Table 3, the global, ensemble mean tropospheric 
ozone column is 30.8 DU, compared to 31.1 DU for 
OMI/MLS, although the latter has a root-mean square 
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Fig. 6. Comparison of the annual mean tropospheric ozone column between (a) the ACCMIP ensemble (different tropopause compared to 
Fig. 2d) and (b) a climatology derived from the Ozone Monitoring Instrument (OMI) and Microwave Limb Sounder (MLS) by Ziemke et 
al. (2011). (c) ACCMIP ensemble bias compared to OMI/MLS (%). See also Tables 2 and 3. 
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Table 3. Tropospheric ozone column, bias and spatial correation for 
the Hist 2000 simulation of the ACCMIP models vs. a climatology 
derived from OMI data (Ziemke et al., 2011). 


Models 

Column/DU 

Bia/DU 

r 

CESM-CAM-superfast 

26.2 

-4.9 

0.79 

CICERO-OsloCTM2 

28.3 

-2.8 

0.84 

CMAM 

30.3 

-0.8 

0.87 

EMAC 

34.8 

3.7 

0.84 

GEOSCCM 

32.1 

1.1 

0.87 

GFDL-AM3 

35.1 

4.0 

0.89 

GISS-E2-R 

33.7 

2.6 

0.85 

GISS-E2-R-TOMAS 

34.7 

3.6 

0.87 

HadGEM2 

28.6 

-2.4 

0.83 

LMDz-OR-INCA 

31.2 

0.1 

0.87 

MIROC-CHEM 

31.3 

0.2 

0.81 

MOCAGE 

28.8 

-2.2 

0.60 

NCAR-CAM3.5 

28.9 

-2.2 

0.84 

STOC-HadAM3 

28.7 

-2.4 

0.81 

UM-CAM 

29.7 

-1.4 

0.75 

ACCMIP mean (± sdev) 

30.8 

-0.4 ±2.7 

0.87 ±0.07 


interannual variability of ~3 DU (Ziemke et al., 2011), which 
overlaps an additional observationally-based estimate from 
the TES instrument of 29.8-32.8 DU (H. Worden, personal 
commnication, 2012). The range from these two instruments 
encompasses the columns calculated by 75 % of the models. 
The spatial correlation between OMI/MLS and the models is 
generally very high (cf. Fig. 6a and b). 

Many of the differences between the ensemble mean 
and OMI/MLS are broadly consistent with the comparison 
against ozonesonde data (Fig. 6c; Table 4), and biases in the 
mean column for a given latitude band are well correlated 
with those for the ozonesondes (r > 0.75 for any pressure 
level). Compared to OMI/MLS, the ensemble mean overes- 
timates the column across the NH mid-latitudes, and under- 
estimates the column over tropical oceans and for all regions 
poleward of approximately 30° S, although the underestimate 
is stronger than suggested by the ozonesonde data. The nega- 
tive bias over the equatorial Pacific in Fig. 6c is not consistent 
with the ozonesonde comparison in Fig. 5, which suggests a 


neutral or positive bias for the mean model. However, this 
region is poorly represented by ozonesonde measurements. 
Correlations between the biases for the NH and SH tropical 
columns are strong (r =0.88), suggesting that similar pro- 
cesses are operating in the regions, even if the sign of the 
bias is different between them (i.e. a model with a stronger 
positive NH tropical bias likely has a SH tropical bias that is 
either positive, or less negative than the ensemble mean). 

Overall, compared to the ensemble of observations, the 
models may have a systematic high bias in the NH, and a sys- 
tematic low bias in the SH. As the emissions are broadly con- 
sistent across the ensemble, the prevalence of this bias could 
suggest they are deficient in some way, in either their amount 
or distribution, or both. However, the models all typically fall 
within the interannual variability of the observations. 

5 Tropospheric ozone from 1850 to 2100 

In this section we consider the changes in tropospheric ozone 
projected by the ACCMIP models for the past (1850 and 
1980), as well as for the near (2030) and more distant (2100) 
future, using the range of RCPs. We begin by discussing 
global-scale changes, followed by regional changes, before 
considering the drivers of the change. 

5.1 Global-scale changes: tropospheric ozone burden 

Figure 7a shows the annual average tropospheric ozone bur- 
den for the ACCMIP models for all the simulations and time 
slices considered. Figure 7b shows the difference in the tro- 
pospheric ozone burden compared to the Hist 2000 simula- 
tion. Data for individual models burdens and their differences 
can be found in Tables 1 and 5 respectively. 

The evolution of the mean tropospheric burden in Fig. 7a 
shows a 25% increase between 1850 and 1980, and a 29% 
increase between 1850 and 2000 (very close to the results 
of Lamarque et al., 2005); the burden increases by 4% be- 
tween 1980 and 2000. Future projections vary with the sce- 
nario. Compared to the mean 2000 burden of 337 Tg, the rel- 
ative changes in the mean burdens for 2030 (2100) for the 
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Table 4. Tropospheric ozone column and bias (DU) for the Hist 2000 simulation of the ACCMIP models vs. the OMI climatology for 
different latitude bands. 



60° S-30° 

S 

30° S-Eq. 

Eq.-30 

° N 

30° N- 

60° N 


Col 

Bias 

Col 

Bias 

Col 

Bias 

Col 

Bias 

CESM-CAM-superfast 

21.9 

-7.8 

21.6 

-8.4 

26.0 

-5.0 

37.2 

3.1 

CICERO-OsloCTM2 

22.1 

-7.6 

27.5 

-2.5 

30.0 

-1.0 

33.2 

-0.9 

CMAM 

29.1 

-0.6 

27.5 

-2.4 

29.3 

-1.7 

36.5 

2.4 

EMAC 

29.4 

-0.3 

33.9 

3.9 

36.7 

5.7 

38.6 

4.5 

GEOSCCM 

28.4 

-1.3 

28.8 

-1.2 

32.0 

1.0 

40.6 

6.5 

GFDL-AM3 

31.4 

1.7 

31.8 

1.9 

35.2 

4.2 

42.9 

8.8 

GISS-E2-R 

27.4 

-2.3 

29.0 

-1.0 

33.5 

2.5 

46.5 

12.4 

GISS-E2-R-TOMAS 

30.0 

0.2 

30.5 

0.5 

33.8 

2.8 

46.3 

12.2 

HadGEM2 

22.8 

-6.9 

26.4 

-3.6 

31.1 

0.1 

34.2 

0.1 

LMDz-OR-INCA 

25.7 

-4.0 

29.2 

-0.8 

31.9 

0.9 

38.3 

4.2 

MIROC-CHEM 

25.1 

-4.6 

31.4 

1.4 

33.8 

2.7 

34.1 

-0.0 

MOCAGE 

18.9 

-10.8 

24.5 

-5.5 

32.3 

1.3 

40.0 

5.9 

NCAR-CAM3.5 

25.1 

-4.7 

24.7 

-5.3 

29.4 

-1.6 

37.6 

3.6 

STOC-HadAM3 

21.5 

-8.3 

26.1 

-3.8 

31.2 

0.2 

35.8 

1.7 

UM-CAM 

26.2 

-3.5 

23.8 

-6.1 

30.4 

-0.6 

40.4 

6.3 

ACCMIP mean (± sdev) 

25.7 ±3.7 

-4.1 

27.8 ± 3.4 

-2.2 

31.8 ± 2.7 

0.8 

38.8 ±4.1 

4.7 

OMI (obs) 

29.7 


30.0 


31.0 


34.1 



Table 5. Differences in the tropospheric ozone burden compared to Hist 2000, using data in Table 1. 


Model 

Hist 


RCP2.6 

RCP4.5 

RCP6.0 

RCP8.5 

1850 

1980 

2030 

2100 

2030 

2100 

2030 

2100 

2030 

2100 

CESM-CAM-superfast 

-111 

-14 

-24 

-72 

- 

- 

-14 

-51 

26 

82 

CICERO-OsloCTM2 a 

-102 

-20 

-11 

-61 

5 

-34 

- 

- 

18 

36 

CMAM 

-84 

-13 

-16 

-57 

7 

-30 

- 

- 

19 

48 

EMAC 

-118 

-25 

- 

- 

1 

-36 

- 

- 

22 

63 

GEOSCCM 

-96 

-13 

- 

- 

- 

- 

- 

- 

- 

- 

GFDL-AM3 

-114 

-8 

-11 

-61 

11 

-22 

12 

-19 

32 

106 

GISS-E2-R 

-92 

-7 

-3 

-40 

9 

-23 

8 

-5 

35 

74 

GISS-E2-R-TOMAS 

-98 

-8 

- 

- 

- 

- 

- 

- 

- 

- 

HadGEM2 

-81 

-18 

-4 

-45 

9 

-12 

- 

- 

23 

70 

LMDz-OR-INCA 

-92 

-17 

-18 

-61 

2 

-29 

-11 

-33 

12 

35 

MIROC-CHEM 

-101 

-20 

-16 

-57 

- 

- 

-3 

-37 

16 

33 

MOCAGE 

-55 

-4 

-3 

-28 

- 

- 

7 

-9 

32 

74 

NCAR-CAM3.5 

-114 

-18 

-19 

-72 

0 

-42 

-14 

-51 

13 

50 

STOC-HadAM3 

-115 

-16 

-19 

-77 

- 

- 

- 

- 

19 

36 

UM-CAM 

-96 

-18 

1 

-28 

16 

1 

- 

- 

29 

75 

Mean 

-98 

-15 

-12 

-55 

7 

-25 

-2 

-29 

23 

60 

Sdev (% of mean) 

17(17%) 

6(39%) 

8(66%) 

16(30%) 

5 (77%) 

13 (52%) 

11 (554%) 

19(65%) 

7 (33 %) 

22 (37%) 


different RCPs are: -4% (-16%) for RCP2.6, 2% (-7%) 
for RCP4.5, -1 % (-9%) for RCP6.0, and 7% (18%) for 
RCP8.5. RCP8.5 is the only scenario to show an ozone 
increase for both time slices (23 Tg and 60 Tg), whereas 
RCP4.5 shows an increase in 2030 (7Tg), before decreas- 
ing in 2100 (— 25 Tg). The ozone burden for the 2030 time 
slice of RCP6.0 is unchanged compared to 2000, although it 
is still higher than 1980. 

Figure 7 also shows a large range in the modelled ozone 
burdens and their differences, with overlapping IQRs be- 


tween many of the time slices. There is a good, but not per- 
fect, correlation (r > 0.7) between the modelled ozone bur- 
den for Hist 2000 and that of other time slices (i.e. models 
generally simulate consistently high or low burdens). How- 
ever, there is no correlation between the modelled ozone bur- 
den and a given burden change, nor between the changes 
in ozone burden for any two periods; i.e. there are no mod- 
els that consistently simulate large (or small) ozone changes 
between time slices, at least at the global scale. This key 
result shows that models are sensitive to emissions and 


Atmos. Chem. Phys., 13, 2063-2090, 2013 


www.atmos-chem-phys.net/13/2063/2013/ 


P. J. Young et al.: Pre-industrial to end 21st century projections of tropospheric ozone 


2077 



1850 1980 2000 2030 2100 2030 2100 2030 2100 2030 2100 
time slice 


■ Historical 

RCP2.6 

12 

RCP4.5 

9 

RCP6.0 

7 

RCP8.5 1 

■v- 

-B 

1 





1850 1980 2030 2100 2030 2100 2030 2100 2030 2100 
time slice 


Fig. 7. (a) Modelled tropospheric ozone burdens for the different 
scenarios and time slices, (b) Change in the tropospheric burden, 
relative to the Hist 2000 simulation. The box, whiskers, line and dot 
show the interquartile range, full range, median and mean burdens 
and differences, and the numbers indicate the number of ACCMIP 
models with results for a given scenario and time slice, all in a sim- 
ilar style to Fig. 1. 


climate changes in different ways. Furthermore, it suggests 
that model weighting schemes based on model skill (e.g. ver- 
sus OMI-MLS) will not necessarily reduce the model spread 
for future projections. A deeper investigation into the drivers 
of this result requires more process-level information from 
the models (e.g. tropospheric ozone budgets from all mod- 
els), and is not possible here. 

The significance of the burden change with respect to Hist 
2000 can also be assessed, using the inter-model spread of 
the differences (Sect. 2.3). This analysis suggests that all 
the changes in the ozone burden are significantly different 
from zero at the 5 % level, except for between Hist 2000 and 
RCP6.0 2030, which is anticipated from Fig. 7b, as this is the 
only time slice where the models do not agree on the sign of 
the change. We again note that “significance” here does not 
mean that the change is significant with respect to interannual 
variability, merely a measure of whether the models agree on 
a change. As shown by Table 5, agreement between models 
on the magnitude of the burden change is generally better for 
the larger changes, namely Hist 1850, RCP2.6 2100, RCP8.5 
2030 and RCP8.5 2100, where, as a percentage of the mean 
change, the standard deviation is 17 %, 30 %, 33 % and 37 % 


respectively. The standard deviations in the differences for 
the other scenarios vary between 40-80 %, although it is very 
large for RCP6.0 2030. 

5.2 Regional-scale changes: burdens, columns and 
concentrations 

Figure 8 shows the mean model regional ozone burden 
changes relative to Hist 2000 for the Hist 1850 and the four 
RCP 2100 simulations, dividing up the troposphere in the 
same manner as Fig. 3b. The figure also indicates the frac- 
tion that each region contributes to the overall ozone change, 
i.e. highlighting asymmetries in the change. From Fig. 7, we 
see that the overall ozone burden change is negative for all of 
these simulations, except RCP8.5 2100. Based on the spread 
of model results, all of the regional burden changes are sig- 
nificantly different from zero. 

For Hist 1850 and RCP2.6 2100 the burden change is 
negative for all regions, with the largest contribution to the 
change coming from the lower ozone precursor emissions in 
the NH extratropics compared to Hist 2000. Unlike for the 
other RCPs, stratospheric ozone recovery (e.g. Eyring et al., 
2010) does not force an increase in tropospheric ozone for 
the SH upper troposphere, despite a 30 % increase in the total 
column ozone (not shown). However, an increase in strato- 
spheric influx is likely masking what would otherwise be 
stronger negative changes due to the precursor decreases (see 
Sect. 5.3). The SH extratropics makes a small contribution to 
the overall change for both the Hist 1850 and RCP2.6 2100 
case. 

The overall decrease for RCP4.5 2100 is about half of 
that between RCP2.6 2100 and Hist 2000, but is still largely 
dominated by the decrease in precursor emissions in the NH 
extratropics, with some contribution from the NH tropical 
lower troposphere. This overall decrease is countered by a 
relatively large increase in the SH upper troposphere, likely 
related to ozone recovery. The magnitude and patterns of ab- 
solute ozone changes are similar for RCP6.0, although the 
tropical upper troposphere makes more of a contribution to 
the overall change than in RCP4.5, in both absolute and rel- 
ative terms. For RCP8.5, ozone increases everywhere, al- 
though the largest contribution is from the 500 to 250 hPa 
pressure band. 

Figures 9, 10 and 11 present information on the annual 
mean spatial patterns of ozone concentration changes, rela- 
tive to Hist 2000, for all the time slices, showing the absolute 
changes in zonal mean ozone, the tropospheric ozone column 
and surface ozone, respectively. Corresponding ozone differ- 
ences for the individual models can be found in the Supple- 
ment (Figs. S4-S6). 

Concentrations for Hist 1850 are less than Hist 2000 ev- 
erywhere except the stratosphere (Fig. 9a), showing the im- 
pact of increased precursors (Fig. 1) and CFC-induced ozone 
depletion respectively. Relative decreases exceed 40 % for 
the column (Fig. 10a) and surface for NH mid-latitudes 
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Fig. 8. ACCMIP ensemble mean change in the tropospheric ozone 
burden relative to the Hist 2000 simulation viewed in Lawrence et 
al. (2001)-style boxes (see also Fig. 3b), for (top to bottom) the 
Hist 1850, RCP2.6 2100, RCP4.5 2100, RCP6.0 2100 and RCP8.5 
2100 simulations. The left hand column shows the absolute dif- 
ference in ozone burden for the different boxes (red/blue for in- 
crease/decrease), with the tropospheric total change indicated in top 
left of each panel. The right hand column shows the fractional (%) 
contribution of each box to the overall change in the tropospheric 
burden. A positive value indicates that the box’s change is the same 
sign as the overall change. Boxes with a fraction larger than 10% 
are highlighted. 


(Fig. 11a). For the latter, absolute decreases are >25ppbv 
for the Mediterranean, much of Asia, and the western USA 
due to less precursor emissions. Differences between Hist 
1980 and Hist 2000 are also significant in many parts of the 
atmosphere. These are distributed in a qualitatively similar 
manner to the Hist 1850 differences, but with smaller de- 
creases due to the smaller change in precursor emissions, al- 
though the lower surface concentrations over South and East 
Asia highlight the recent emission growth in that region (e.g. 
Zhang et al., 2009). A notable non-significant change is seen 
for surface ozone concentrations over the eastern US, where 


~50% of the models simulate higher ozone for 1980. This 
is in qualitative agreement with the recent analysis of ozone 
trends by Parrish et al. (2012), although transient simula- 
tions, better constrained to observed interannual changes in 
meteorology and emissions would be needed to investigate 
this further. 

For RCP2.6, the distribution of tropospheric ozone 
changes is very similar to the Hist 1850 difference, with peak 
reductions of 30-40% by 2100 (Figs. 9c, 10c and 11c) re- 
flecting the partial reversal of the anthropogenic ozone pre- 
cursor changes compared to 1850 and 2000. Similar pat- 
terns are evident again for RCP4.5 and RCP6.0 in 2100 
(Fig. 9f and h), although (non-significant) increases in SH 
ozone penetrate deeper into the troposphere compared to 
RCP2.6 (see below). RCP4.5 also shows significant increases 
throughout the tropical and SH troposphere for the 2030 time 
slice (Fig. 9e), reflecting a redistribution in precursor emis- 
sions: tropical NO x emissions are higher for this scenario in 
2030 than all others, except RCP8.5. Despite non-significant 
changes in the ozone burden for RCP6.0 2030 (Fig. 7), there 
are significant decreases in ozone in the tropical regions 
(Fig. 9f). RCP8.5 has significant increases throughout the 
troposphere, except for surface concentrations and the col- 
umn over the equatorial Pacific in 2100, where the dominant 
effect may be increased specific humidity in the warmer cli- 
mate increasing the ozone loss rate (e.g. Johnson et al., 1999; 
Zeng and Pyle, 2003). 

The zonal mean ozone changes by 2100 for the different 
RCPs are qualitatively similar to those presented by Kawase 
et al. (2011), except that the ACCMIP ensemble mean does 
not show upper tropospheric ozone increases in the NH mid- 
latitudes for RCP4.5 and 6.0. Sensitivity experiments by 
Kawase et al. (2011) demonstrated the importance of an en- 
hanced Brewer-Dobson circulation (BDC) and recovery of 
stratospheric ozone levels in increasing future upper tropo- 
spheric ozone levels for RCP4.5, 6.0 and 8.5 in 2100. Strato- 
spheric chemistry-climate models are robust in projecting re- 
covery of the ozone layer due to a reduction in halogen levels 
(Eyring et al., 2010), as well as an intensification of the BDC 
with increased greenhouse gas concentrations (Butchart et 
al., 2006). The change in zonal mean ozone in Fig. 9 is char- 
acteristic of these processes, and, in particular, the reduc- 
tion in tropical lower stratospheric ozone and enhancement 
of high latitude ozone is indicative of stronger BDC (Randel 
et al., 2002). This tropical/li igli-1 ati tudc seesaw pattern of the 
changes intensifies with the increased climate change from 
RCP2.6 to 8.5, further illustrating the coupling of the BDC 
to greenhouse gas levels in these models (see also Lamar- 
que et al., 2011). The mid-latitude peaks in the tropospheric 
ozone column changes (Fig. lOj) are also indicative of in- 
creases in the stratospheric ozone influx (Olsen et al., 2004). 
For RCP8.5, the simulations of Kawase et al. (2011) also 
showed that the very large increase in methane levels was 
driving increases throughout most of the troposphere (see 
also Brasseur et al., 2006; Fiore et al., 2008). 
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Fig. 9. Absolute change in annual zonal mean ozone for the AC- 
CMIP ensemble mean compared to the Hist 2000 simulation (ppbv). 
Top row shows the difference for the Hist 1850 and 1980 time slices. 
The next four rows show the difference for the 2030 and 2100 time 
slices of the RCP simulations. Non-white regions indicate where 
the change is significant at the 5 % level, based on the spread of the 
differences between the models. The red dashed line indicates the 
position of the annual zonal mean 150 ppbv ozone contour from the 
Hist 1850 simulation. 


With the exception of the upper troposphere and lower 
stratosphere, the ozone differences in regions with non- 
significant changes are small, and it may be that the mod- 
els agree that these changes in these regions are not signif- 
icant in the context of interannual variability (Tebaldi et al.. 
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Fig. 10. As Fig. 9, but for the absolute change in the tropospheric 
ozone column (DU). 


2011). For the regions of significant change, the standard de- 
viation of the differences can exceed 100% of the ensem- 
ble mean difference, though generally only for surface and 
column values close to emission regions. (By construction, 
the colour-filled regions of Figs. 9-11 indicate where most 
models agree on the sign of change, but they do not show 
where there are large ranges modelled for positive or neg- 
ative changes.) Flowever, as noted for changes in the total 
tropospheric burden, there is no apparent correlation with a 
model’s present day ozone level in one region and the change 
in ozone that is modelled for the same region. 
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Fig. 11. As Fig. 9, but for the absolute change in surface ozone 
(ppbv). 


Model agreement on the distribution of the differences 
is very good, with most models being highly spatially cor- 
related with the ensemble mean difference (see Figs. S4- 
S6). Notable exceptions are LMDz-OR-INCA and CICERO- 
OsloCTM2, which have fixed stratospheric ozone levels or 
influx. This testifies to the potential importance of strato- 
spheric circulation and ozone changes for tropospheric ozone 
projections, as most of those models that include some repre- 
sentation of stratospheric ozone evolution have changes that 
are generally well correlated with each other. Flowever, for 
RCP8.5 2100 MOCAGE and STOC-HadAM3 are not well 
correlated with the other models, despite having some repre- 


sentation of stratospheric ozone change. MOCAGE strongly 
concentrates tropospheric ozone increases at high latitudes 
rather than the mid latitudes and tropics as per Figs. 9j and 
lOj. This could relate to systematic differences in the lo- 
cations of stratosphere-to-troposphere ozone transport, the 
distribution of stratospheric ozone, or a combination of the 
two. While STOC-FIadAM3 makes use of the time-evolving 
stratospheric ozone dataset of Cionni et al. (2011), it uses 
it to help constrain ozone concentrations at the model top, 
rather than simply to overwrite ozone above the tropopause. 
This may account for the outlying stratospheric ozone trends 
seen for this model in the supplementary material. A STOC- 
FIadAM3 simulation where only the climate is changed to 
RCP8.5 2100 conditions concentrates tropospheric ozone 
column increases over continental regions, perhaps implying 
a stronger role for LNO x increases instead of more strato- 
spheric ozone influx (see supplementary material of Steven- 
son et al. (2012) and below). 

5.3 Ozone budget changes for a subset of models 

We can gain additional insight into some of the processes 
controlling global tropospheric ozone changes from the lim- 
ited ozone budget data. Figure 12 shows the percentage 
change in the tropospheric ozone budget terms (P, L, P mi- 
nus L, D, S; n f and r) for all the scenarios and time slices 
compared to Hist 2000, showing only the five ACCMIP mod- 
els where there are sufficient data (see Table 2; GEOSCCM 
had only Hist data available). We concentrate on the relative 
changes to minimise the impact of different budget defini- 
tions and the range of different model scheme complexities 
(e.g. higher VOC emissions tend to mean higher P and L). 
Note, not all five models are represented for each variable or 
time slice. 

Figure 12a-b show that the relative changes in the P and L 
terms qualitatively resemble the changes in the tropospheric 
ozone burden in Fig. 7b. Individual models tend to agree on 
the magnitude of the relative changes to within 10-20%, 
although absolute changes differ more. Compared to Hist 
2000, changes in the net chemical production (NCP, P minus 
L; Fig. 12c) are overwhelmingly negative for all time slices 
and models (the RCP8.5 2030 time slice change for STOC- 
HadAM3 is the only exception), likely aided by an increase 
in the water-mediated loss of ozone (via 0( 1 D) + H 20 ) for 
the RCPs due to higher specific humidity in the warmer cli- 
mates (see Fig. 10 of Lamarque et al., 2013). Despite the 
reductions in NCP, for most models the absolute value of 
NCP is positive for all time slices, with net chemical destruc- 
tion only shown for some CESM-CAM-superfast (Hist 1850, 
RCP2.6. 2100 and RCP6.5 2100) and UM-CAM (Hist 1850, 
RCP2.6. 2100, RCP4.5 2100 and RCP8.5 2100) time slices. 
The relative changes in D (Fig. 12d) are qualitatively similar 
to those for P and L, although, notably, with smaller relative 
changes for RCP8.5 2100. Changes in D will depend on the 


Atmos. Chem. Phys., 13, 2063-2090, 2013 


www.atmos-chem-phys.net/13/2063/2013/ 






P. J. Young et al.: Pre-industrial to end 21st century projections of tropospheric ozone 


2081 


(a) AP 


(d) AD 


■ Hist 

RCP2.6' 

RCP4.5 

RCP6.o' 

RCP8.SQ ■ 

1 1 

— i — i — i — i — i — i — i — i — i — i — ■ — 

- two 

o 

§ 

8 

8 

oo< 

* 

a 

o 

a 



Hist 

RCP2.6 

RCP4.5 

A 

RCP6.0 

RCP8.5 

g : 

o ■ 

8 

0 

ft 

A 

A O 

o 

6 o 



° & 

Q 

A 

■ 


Q 


o 


■ ■ ■ i 

D»» 

O 





1850198020302100203021002030 2100 20302100 
time slice 


1 850 1 980 2030 21 00 2030 21 00 2030 21 00 2030 21 00 
time slice 



Hist 

RCP2.6' 

RCP4.5 

o . 

RCP6.0 

RCP8.5 

0 

0 

o 0 

o 

o 

o 

- 

6 

9 . 

o 





1 850 1 980 2030 21 00 2030 21 00 2030 21 00 2030 21 00 
time slice 


1 850 1 980 2030 21 00 2030 21 00 2030 21 00 2030 21 00 
time slice 


(C) ANCP 

40 
20 
0 

-20 

1 850 1 980 2030 21 00 2030 21 00 2030 21 00 2030 21 00 
time slice 

|®CESM-CAM-superfast ^GFDL-AM3 ANCAR-CAM3.5 ^STOC-HadAM3 ©UM-CAM 



(f)At 



time slice 


Fig. 12. Relative change in tropospheric ozone budget terms compared to Hist 2000, for all scenarios and time slices, for the subset ACCMIP 
models with data. Construction is similar to Fig. 7b, except the box/whisker is replaced with a symbol for each model. Relative changes are 
shown for (a) chemical production (P), (b) chemical loss (L), (c) net chemical production (NCP; P minus L), (d) dry deposition (D), (e) 
inferred stratospheric input (S; n f) and (1) lifetime (r). 


distribution of near-surface ozone changes (Fig. 11) as well 
as the characteristics of the deposition schemes. 

Figure 12e shows that changes in S m r vary more between 
the models. Moreover, the size of the change is qualitatively 
related to the magnitude of S; n f in the Hist 2000 slice: UM- 
CAM has the largest Hist 2000 S m f and the largest changes in 
Sinf, whereas STOC-FladAM3 has the smallest Hist 2000 S; n f 
and the smallest changes in S; n f (see Sect. 5.2 for further dis- 
cussion on the different treatment of the stratosphere). Except 
for some STOC-HadAM3 results (Hist 1850 and RCP2.6 
2100), all models show an increase in Sj n f for all time slices 
compared to Hist 2000. As most simulations also show a de- 
crease in P, the increases in Si n f point to the increased im- 
portance for that tropospheric ozone source term under his- 
torical and projected conditions. For UM-CAM, S; n f is half 
the size of P for Hist 1850, 40 % of P for RCP2.6 2100 and 
30% of P for RCP8.5 2100; S; n f varies between 7-22% of 


P for the other models. The general increase in S; n f for the 
RCPs is consistent with the qualitative analysis in Sect. 5.2 
and Kawase et al. (2011). 

Figure 12f shows that the relative changes in r are smaller 
than for the other terms, being most notable for the scenar- 
ios with larger emission reductions compared to Hist 2000. 
These are namely Hist 1850 (3.6 to 7.8 day longer lifetime) 
and, to a lesser extent, RCP2.6 2100 (0.2 day shorter life- 
time to 4.7 day longer lifetime). While Kawase et al. (2011) 
are mostly in agreement with this limited ACCMIP ensemble 
for RCP2.6. and RCP8.5 lifetime changes, they report life- 
time decreases for RCP4.5 and RCP6.0, where the ACCMIP 
model spread suggests an ambiguous result. 

Overall, as in Sect. 4.1, we caution that the ACCMIP re- 
sults cannot be used for a consistent comparison, due to the 
different methods of determining P and L (and therefore 


www.atmos-chem-phys.net/13/2063/2013/ 


Atmos. Chem. Phys., 13, 2063-2090, 2013 


2082 


P. J. Young et al.: Pre-industrial to end 21st century projections of tropospheric ozone 


(b) 

- Hist 
RCP2.6 
RCP4.5 

RCP6.0 

RCP8.5 

H 

H 

1 

r 



-40 -30 -20 -10 0 10 -2000 0 2000 4000 


ATotal NO x (Tg N/yr) 


ACH 4 (Tg) 


Fig. 13. ACCMIP ensemble mean change in the tropospheric ozone burden compared to the Hist 2000 simulation as a function of (a) 
changes in total NO x emissions and (b) changes in the tropospheric methane burden. Error bars indicate ± 1 std. dev. of the changes in 
ozone, NO x emissions and methane burdens, calculated from the spread of the models. Different colours represent the different scenarios, 
whereas different symbols represent the different time slices. 


NCP), blit are nevertheless instructive for future compar- 
isons. 

5.4 NO x , methane and the implied role of climate 
change as drivers of the total ozone changes 

The relationship of modelled tropospheric ozone burdens 
with methane and NO x is well established (Stevenson et 
al., 2006; Wild, 2007; Fiore et ah, 2008; Wild et ah, 2012) 
and Fig. 13 presents how the ACCMIP ensemble mean 
ozone burden changes for each simulation, together with (a) 
changes in the mean NO x emission and (b) changes in the 
mean tropospheric methane burden. 

Figure 13a shows that the evolution of tropospheric ozone 
from the pre-industrial period to present day tracks the 
change in NO x emissions in a near linear fashion, similar 
to the relationship presented by Stevenson et ah (2006). The 
decrease in NO x emissions for R.CP2.6, 4.5 and 6.0 sees the 
ozone burden decrease again, although at a slightly reduced 
rate than for the Hist simulations partially due to the redistri- 
bution of precursor emissions equatorward, where the ozone 
production efficiency is greater (Gupta et ah, 1998; Wild and 
Palmer, 2008). This is particularly the case for RCP4.5 2030, 
which - as noted in Sect. 5.2 - has an increase in tropical 
NO x emissions compared to Hist 2000, despite an overall de- 
crease. 

RCP8.5 is the clear outlier for the simple NO x -ozone re- 
lationship, with a 60 Tg increase in the tropospheric ozone 
burden coupled with a 12 TgN yr -1 reduction in NO x emis- 
sions. Based on the results of the other simulations, and only 
considering NO x changes, we might expect a 40-50 Tg de- 
crease in the tropospheric ozone burden. However, as already 
stated, a defining feature of RCP8.5 is the large increase in 


methane concentrations through the 21st century, and the re- 
lationship between ozone changes and methane changes for 
the simulations is shown in Fig. 13b. Taking in data across all 
the simulations shows that the relationship is not linear, and it 
clearly partially depends on the levels of other ozone precur- 
sors, as well as their distribution (e.g. see Wild, 2007). For 
instance, the change between Hist 2000 and RCP8.5 2030 
qualitatively sits on the same line as the ozone-methane re- 
lationship for the Hist simulations, likely due to the methane 
increase in RCP8.5 2030 being accompanied by an increase 
in NO x emissions. Between 2030 and 2100, the reduction 
NO x (and other) emissions for RCP8.5 contributes to the fact 
that a given methane increase does not produce as much of 
an ozone increase. 

The impacts of climate change further complicate this cor- 
relation. Using the parameterised relationship between ozone 
abundance and the levels of its precursor emissions (but ex- 
cluding climate) developed by Wild et al. (2012), we would 
expect an ozone burden increase of approximately 30 Tg be- 
tween Hist 2000 and RCP8.5 in 2100. However, at 60 Tg the 
ACCMIP ensemble mean increase in ozone is double that 
expected, and consistent with equal roles for methane in- 
creases and the net impacts of climate change, i.e. through 
promoting increased influx of stratospheric ozone, chang- 
ing LNO x , and impacting reaction rates, through tempera- 
ture and humidity changes. This is broadly similar to re- 
sult from the RCP8.5 sensitivity simulations of Kawase et 
al. (2011), which showed a 5.5 DU increase in the global 
mean tropospheric column when all drivers changed, and a 
2.0 DU increase when only (non-methane) greenhouse gases 
changed - i.e. assuming linearity, climate change accounted 
for 36% of the total tropospheric ozone change. We have 


Atmos. Chem. Phys., 13, 2063-2090, 2013 


www.atmos-chem-phys.net/13/2063/2013/ 


P. J. Young et al.: Pre-industrial to end 21st century projections of tropospheric ozone 


2083 


some evidence for the importance of the stratosphere in 
the ACCMIP ensemble from the large increases in S m r for 
RCP8.5 2100 (Fig. 12d), although it would be more instruc- 
tive to have this data for all models. 


6 Discussion 

Considering the full ACCMIP ensemble, the results for 
ozone change are most unambiguous for Hist 1850, RCP2.6 
2100 and RCP8.5 2100, both in terms of magnitude and 
distribution. These represent the extremes of the spectrum 
of historical and RCP scenarios, with the former two hav- 
ing the lowest concentrations of all ozone precursors, and 
RCP8.5 having relatively low NO x , CO and NMVOC emis- 
sions, but very high methane coupled with a strong warming 
(see Lamarque et ah, 2013). With the generally low concen- 
trations of the more “complex” VOC precursors, “basic” tro- 
pospheric chemistry (i.e. involving methane, CO, NO x , HO x 
and ozone) becomes more important. The reactions describ- 
ing this chemistry are generally very similarly represented 
in most models (e.g. Emmerson and Evans, 2009), and, with 
the reduced importance of the chemistry of more complex 
VOCs, this could potentially be driving a lot of the similar- 
ity between the models. The relative changes in the P and 
L tenus for the subset of ACCMIP models with budget data 
are reasonably well clustered (Figs. 12a and 12b), but un- 
certainty in the interpretation would be reduced if the whole 
ensemble were better represented. Further useful information 
could come from a systematic investigation of the response 
of ozone to idealised precursor changes in the different mod- 
els, such as through the sensitivity studies of Wild (2007). 

While there is good agreement between the models for the 
ozone changes between Hist 1850 and Hist 2000, we note 
that none of the ACCMIP models can reproduce the low 
surface ozone concentrations suggested by late- 19th century 
measurements using the Schonbein method (Pavelin et ah, 
1999; Hauglustaine and Brasseur, 2001). Compared to the 
data presented by Hauglustaine and Brasseur (2001), biases 
for the ensemble mean are 40-350% (not shown). This re- 
sult has not changed greatly over the last two decades (e.g. 
Pavelin et ah, 1999). Indeed, one of the only model studies 
to simulate ozone in-line with these Schonbein data invoked 
large perturbations in the emissions of VOC and NO x , com- 
pared to those imposed in this study (Mickley et ah, 2001). 

As well as uncertainty in ozone precursor emissions, there 
is scope for uncertainty in the representation of the oxidation 
chemistry during the cleaner pre-industrial period, where, in 
particular, levels of NO x are expected to have been much 
lower than today. Since isoprene emissions may not have 
changed dramatically since the pre-industrial period (e.g. Ar- 
neth et ah, 2010), modifications to the low-NO x chemistry of 
isoprene may be important. Novel isoprene chemistry has re- 
cently been included in simulations of pre-industrial ozone 
by Archibald et ah (2011). However, the changes they im- 


posed to rectify problems with simulating surface OH in the 
tropics led to an increase in surface ozone everywhere. Bio- 
genic hydrocarbon emissions play a dual role in ozone pro- 
duction acting as ozone precursors on the one hand but also 
many are able to react directly with ozone at a fast rate, or 
decrease ozone production by sequestering NO x (Fiore et ah, 
2005; Horowitz et ah, 2007; Young et ah, 2009). Recently, 
tropospheric halogen chemistry has been postulated as be- 
ing an important process missing in many models that have 
attempted to simulate pre-industrial ozone (Parrella et ah, 
2012; Saiz-Lopez et ah, 2012), but such processes were not 
included in the ACCMIP models. Clearly, more understand- 
ing of all these processes is important for simulating past and 
future tropospheric ozone. 


7 Summary and conclusions 

This study has analysed tropospheric ozone changes from 
1850 to 2100 from the range of chemistry models that con- 
tributed to the Atmospheric Chemistry and Climate Model 
Intercomparison Project (ACCMIP), running the latest set 
of ozone precursor emissions scenarios, and with 14 out of 
15 models also including representations of the changing 
climate. The ensemble mean ozone distribution compares 
favourably with present day satellite and ozonesonde obser- 
vations. The seasonal cycle is well captured, except in some 
locations in the tropical upper troposphere, and there are sug- 
gestions of a high bias in the Northern Hemisphere (NH) and 
a low bias in the Southern Hemisphere (SH), although most 
model results fall within the range of observed interannual 
variability. 

In agreement with other studies (e.g. Lamarque et ah, 
2005), the modelled tropospheric ozone burden in 1850 is 
~30 % lower than the present day, with the largest contribu- 
tion to the change coming from the NH extratropics. Inter- 
model agreement on the magnitude of this change is rea- 
sonably high (98 ± 17 Tg), although modelled surface ozone 
concentrations are higher than the available pre-industrial 
measurements (as per Hauglustaine and Brasseur, 2001) sug- 
gesting that there are still unresolved issues with correctly 
modelling pre-industrial ozone levels (Mickley et ah, 2001). 
Modelled ozone also increases somewhat between 1980 and 
2000, particularly over industrialised regions in the NH, in 
agreement with the general picture described by Parrish et 
ah (2012). 

Future changes in tropospheric ozone were considered for 
2030 and 2100 time slices, using projections of climate and 
ozone precursor emissions from four Representative Con- 
centration Pathways (RCPs). Compared to 2000, the relative 
changes for the tropospheric ozone burden in 2030 (2100) 
for the different RCPs are: -4% (-16%) for RCP2.6, 2% 
(-7%) for RCP4.5, -1% (-9%) for RCP6.0, and 7% 
(18%) for RCP8.5. The decreases apparent for most RCPs 
are due to reductions in precursor emissions, but the increase 
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in ozone for RCP8.5 is in spite of reductions in nitrogen 
oxide emissions. From the limited ACCMIP ozone budget 
data and as implied by comparison to other studies (Kawase 
et ah, 2011; Wild et ah, 2012), the RCP8.5 ozone increase 
can be attributed to the very large increase (~ doubling) in 
methane and increased stratospheric influx of ozone (40- 
150% increase, as determined from the five models with 
ozone budget data). Inter-model agreement on the magni- 
tude of the total change compared to 2000 is best where the 
changes are large, such as for RCP 2.6 in 2100 (—55 ± 16 Tg) 
and RCP8.5 in 2100 (60 ± 22 Tg). While models with higher 
present day ozone burdens have higher ozone burdens for 
the other time slices, there is no relationship between bur- 
den and burden change, or between burden changes for dif- 
ferent periods/scenarios. This key result suggests that any 
model weighting schemes based on the present day model 
bias (e.g. versus OMI-MLS) will not necessarily reduce the 
model spread for future projections. 

For the changes with all RCPs, generally the models are 
highly spatially correlated with one another, agreeing on 
where the changes are occurring, if not their magnitude. No- 
table exceptions are for the models that do not include repre- 
sentation of the changing influence of the stratosphere on the 
troposphere, highlighting the importance of ozone recovery 
(Eyring et ah, 2010) and climate change-induced circulation 
strengthening (Butchart et ah, 2006; Butchart et ah, 2010; 
SPARC-CCMVal, 2010) in this region of the atmosphere. 

Overall, we have shown that the multi-model mean of AC- 
CMIP models generally simulates present day tropospheric 
ozone well, and agrees on the sign of past and future tro- 
pospheric ozone changes and how those changes are dis- 
tributed. Fully establishing the consistency of the chemi- 
cal and physical processes driving the changes in the dif- 
ferent models is limited in this study due to lack of com- 
parable ozone budget statistics. For example, we cannot ex- 
plain the lack of a correlation between the present day tropo- 
spheric ozone burden and pre-industrial to present day bur- 
den change, something which may be attributed to different 
model sensitivities to precursor emission changes, combined 
with a range of ozone fluxes from the stratosphere. Future 
studies will require careful thought as to how to make the 
necessary diagnostics comparable across models and we en- 
courage discussion and further study within the chemistry- 
climate community of how best to manage this. In addition, 
the chemistry community would benefit from a thorough in- 
vestigation into the importance of the different processes that 
control tropospheric ozone, as has been examined for cloud 
condensation nuclei in aerosol models (Lee et ah, 2012). 
Finally, like Kawase et al. (2011), this study highlights the 
strong influence of stratospheric processes in controlling tro- 
pospheric ozone, which may encourage more tropospheric 
chemistry-climate modelling groups to move towards full at- 
mosphere simulations. 
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